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ABSTRACT 

We present a theoretical and exact analysis of the bispectrum of projected galaxy cata- 
logues. The result can be generalized to evaluate the projection in spherical harmonics 
of any 3D bispectrum and therefore has applications to cosmic microwave background 
and gravitational lensing studies. 

By expanding the 2D distribution of galaxies on the sky in spherical harmonics, we 
show how the 3-point function of the coefficients can be used in principle to determine 
the bias parameter of the galaxy sample. If this can be achieved, it would allow a 
lifting of the degeneracy between the bias and the matter density parameter of the 
Universe which occurs in linear analysis of 3D galaxy catalogues. In previous papers we 
have shown how a similar analysis can be done in three dimensions, and we show here 
through an error analysis and by implementing the method on a simulated projected 
catalogue that ongoing three-dimensional galaxy redshift surveys (even with all the 
additional uncertainties introduced by partial sky coverage, redshift-space distortions 
and smaller numbers) will do far better than all-sky projected catalogues with similar 
selection function. 

Key words: cosmology: theory - galaxies: clustering - bias - large-scale structure of 
the Universe 



1 INTRODUCTION 

The clustering of mass in the Universe is an important fossil record of the early perturbations which gave rise to large-scale 
structure today. Knowledge of the mass clustering puts powerful constraints on the quantity and properties of dark matter in 
the Universe, and the generation mechanism for the perturbations. Most of our knowledge of the mass clustering is, however, 
indirect, coming principally from the distribution of galaxies. A major obstacle in interpretation is therefore the uncertain 
relationship between galaxy and mass clustering - a relationship which is conventionally quantified by the 'bias parameter', 
6. In particular, attempts based on linear perturbation theory to measure the density parameter of the Universe, Qq, through 
peculiar velocity or redshift-distortion studies, yield only the degenerate combination (3 = Q,®' 6 /b, making it impossible to 
determine Oo without determining b. The degeneracy can be lifted by going to second order in perturbation theory, and this 
can be achieved most elegantly by studying the bispectrum, which is the three-point function in Fourier (or spherical harmonic) 
space. A major positive feature of the bispectrum method is that it can provide error bars on the desired parameters. The 
method works because gravitational instability leads to a density field which is progressively more skewed to high densities 
as it develops, and this skewness appears as a non-zero bispectrum. This behaviour can also be mimicked by biasing, if the 
galaxy density field is a local, nonlinear function of the underlying mass density field. This possibility must be dealt with. 
The two effects can, however, be separated by the use of information about the shape of the structures: in essence the effect 
of biasing is to shift iso-density contours up (or down), while maintaining the shape of the contour; gravitational evolution, 
instead, changes the shape, usually leading to flattening of collapsing structures (e.g. Zel'dovich pancakes). 

Fry (1994) recognized the role the bispectrum could play in determining the bias parameter, and Matarrese, Verde 
& Heavens (1997) (hereafter MVH97) and Verde et al. (1998) have turned the idea into a practical proposition for 3D 
galaxy redshift surveys by including analysis of selection functions, shot noise, and redshift distortions (see also Scoccimarro, 
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Couchman & Frieman 1999). The latter is potentially a serious problem for 3D surveys, as the signal for bias comes from mildly 
nonlinear scales, where the redshift distortions are not trivial to analyze. However, experiments on simulated catalogues (Verde 
et al. 1998) show that the method is successful. Note however that the theory has been developed only in the 'distant-observer 
approximation' (see e.g. Kaiser 1992), and is applicable to relatively deep surveys such as the Anglo- Australian two-degree 
field galaxy redshift survey (Colless 1996; Colless 1999) and the Sloan Digital Sky Survey (Gunn 1995). Shallow surveys such 
as the IRAS PSCz should be analysed in spherical coordinates (cf Taylor & Heavens 1995; Tadros et al. 1999 for the power 
spectrum), and suffer from high shot noise, so cannot be usefully used for bispectrum analysis based on a Fourier expansion 
(see also Bharadwaj 1999). 

The absence of suitable existing 3D surveys prompts us to consider whether the bias might be extracted from a projected 
galaxy catalogue, such as the APM galaxy survey (Maddox et al. 1990; Loveday et al. 1992; Loveday 1996). With only angular 
positions, the information is more limited, but the survey is not complicated by redshift distortions, and contains a large 
number ~ 10 6 of galaxies. The DPOSS catalogue (Djorgovski et al. 1998) will be even larger, with 50 million galaxies and 
nearly all-sky. There are two important caveats to keep in mind: first, to have a measurement of Qo we need to be able to 
measure the (5 parameter and the linear bias parameter. It is not possible to extract the f3 parameter from a two-dimensional 
survey, /3 will need to be determined from a -different- three-dimensional survey. The selection criteria will necessarily be 
different for different catalogues, and so will be the galaxy population selected. Since different galaxy populations can have 
different bias with respect to the underlying dark matter distribution, some care needs to be taken in the interpretation of 
the final result i.e. the value for Q.®. The other caveat concerns the effect of the evolution along the line of sight (also referred 
to as the light-cone effect). This is due to the fact that galaxy clustering evolves gravitationally with time along the line of 
sight and depends on the (unknown) cosmology. For shallow surveys such as the APM, this effect is smaller or comparable to 
the cosmic variance, in what follows we will neglect this effect for this reason. However, for deeper surveys, this effect needs 
to be properly taken into account. Assuming these issues can be dealt with, the key requirement is to obtain an expression 
for the projected bispectrum given an analytical formula for the spatial one. An expression for the projected bispectrum in 
the small angle approximation has been presented by Buchalter, Kamionkowski & Jaffe (1999). However this might not be a 
good approximation for the bispectrum if the sample is close enough to the observer or if the scales under analysis are large. 
In fact it is not known a priori whether, for the bispectrum, the small angle approximation is valid on large enough scales for 
the second order perturbation theory to hold: it is necessary to obtain an exact expression for the projected bispectrum using 
spherical harmonics expansion. Only then it will be possible to test the limit of validity for the small-angle approximation 
bispectrum (Verde et al. 2000). 

In this paper, we develop the theory for projected catalogues in a full treatment. The resulting expression for the spherical 
harmonic projected bispectrum can straightforwardly be applied to gravitational lensing and to cosmic microwave background 
(CMB) studies, for comparison with observations such as the claimed detection by Ferreira, Magueijo & Gorski (1998). 

In Section 2 we expand the sky density of galaxies in spherical harmonics with coefficients a™, and compute an explicit 
expression for the bispectrum (a^a^a™ 3 ) accurate to second-order in perturbation theory. In particular, we show how this 
quantity depends on the bias parameter. We present in Section 3 an error analysis specific to the second-order perturbation 
theory bispectrum, which shows the expected uncertainty in the derived bias parameter, and test on a numerical simulation. In 
the Appendices, we detail asymptotic results which are useful for high-i? spherical harmonics. The main conclusion of large-scale 
structure application in this paper is that 3D large scale structure surveys (even with small sky coverage, smaller numbers, 
and the complications of redshift-space distortion, shot noise etc.) will do far better than all-sky projected catalogues for the 
purpose of measuring the bias parameter. However the mathematics developed for this purpose has much wider applications: 
with appropriate radial weight functions, the analysis can be applied to the CMB bispectrum induced by lensing, Sunyaev- 
Zel'dovich effect, the integrated Sachs- Wolfe effect or foreground point sources, and to gravitational lenses studies. 

2 PROJECTED BISPECTRUM IN SPHERICAL HARMONICS 

Let the projected galaxy density field be n(f2), where f2 represents angular positions in the sky. If the three-dimensional galaxy 
density field is p(r) (with mean p) and the selection function is tp(r), the projected density is 




(1) 



We expand the projected density in spherical harmonics (see Appendix A for definitions) 




(2) 



where 8 = (p — p)/p is the fractional overdensity in galaxies. The average surface density is 
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n = J dr r 2 ptp(r). (3) 

and is inserted in the transform for convenience. The three-point function of the coefficients may be factorised by isotropy 
(e.g. Luo 1994) into 

/ mi m-2 7713 \ D /^l ^2 ^3 \ ( A\ 

\ Cl l 1 a t 2 a t 3 I — \m im2 m 3 ) (4J 

where (m^ 2 ,^ 3 ) is the Wigner 3J symbol. We refer to Bi 1 i 2 £ 3 as the angular bispectrum. From general considerations about 
rotational invariance of the quantity {a^a^a^) the indices £i,rm for i = 1,2,3 must satisfy the following conditions: 

(i) £j +£ k >li>\tj-lk\ (triangle rule) 
(ii) l\ + £2 + £3 — even 
(in) mi + Wi2 + rn,3 = 0. 

The presence of the 3J symbol ensures that these conditions are satisfied. 

In order to be able to extract the bias parameter from projected catalogues, the effect of the projection in the configuration 
dependence of the bispectrum needs to be understood (Fry & Thomas 1999) . 

To do so, we compute the angular bispectrum in terms of the 3D bispectrum, B(ki,k 2 ,k 3 ) defined by (<5ki<5k2<5k3) = 
(2-7r) 3 _B(ki, k 2 , k 3 )5 I5 (ki + k 2 + k 3 ), where the Fourier transform of 5 is 5k = J d 3 r <5(r) exp(ik • r), and 5 D denotes the Dirac 
delta function. 

We proceed from (Q): 

(a™ 1 ^ 2 ^ 3 ) = (S) 3 / ^i^2dn 3 dr 1 dr2dr 3 r^|rf^ 1 V'2V'3(5(r 1 )5(r 2 )5(r3))Y^ 1 (fi 1 )Y«7(n 2 )y^ 3 (n3) . (5) 
The 3D three-point function (in real space) is related to the 3D bispectrum by 

(*(ri)a(pa)«(rs)) = [ d 3 k 1 d 3 k 2 d 3 k 3 B(k 1 , k 2 , k 3 )e l(kl • ri+k2 ' r2+k3 ' r3) 5 D (ki + k 2 + k 3 ) . (6) 



(7) 



(2tt) 6 j 

We then define the quantity: 

J(n,r 2 ,r 3 )= / dkidk 2 dk 3 kiklk 2 dtt kl dSl k3 dn ka B(ki, k 2 , k 3 )e l(kl ' ri+k2 r2+ka r3) 5 D (ki + k 2 + k 3 ) 

JO J 4vr 

because we will later expand the exponential in spherical harmonics and perform the angular integrations in (Q) explicitly. 
In second order perturbation theory the bispectrum is: 

B(ki,k 2 ,k 3 ) = /C(ki,k2)P(fei)P(fe)+c2/c. , (8) 

where the shape-dependent factors K, can be found in Fry (1984), Catelan et al. (1995) and MVH97. The dependence of K. 
on the cosmology in negligible (e.g. Kamionkowski & Buchalter (1999) and references therein), so in what follows we assume 
an Einstein-de Sitter Universe. The factors are, however, dependent on the biasing model assumed. If we take a local biasing 
model, then for consistency with second-order perturbation theory, we expand in a Taylor series the galaxy overdensity to 
second-order in the matter overdensity <5 m : 

S(x) = M m (x) + i& 2 «5 m (x) 2 , (9) 

(a constant term bo is irrelevant except at k = and is ignored) . Here 61 is the linear bias parameter and & 2 is the quadratic 
bias parameter. The linear bias parameter b that appears in the definition of f3 and that is needed to recover Qo, is b = 61 
on large scales, under fairly general conditions (Heavens, Matarrese & Verde 1998). Note that we take the bias function to 
be deterministic, not stochastic (cf Cen & Ostriker 1992; Dekel & Lahav 1999; Tegmark & Bromley 1999; Matsubara 1999; 
Somerville et al. 1999); it has been shown (Taruya, Koyama & Soda 1999) that the effect of stochastic bias on the bispectrum 
is very similar to that of nonlinear bias (Eq. ^). This formalism might be straightforwardly extended to the case when the 
bias process operates in Lagrangian, rather than Eulerian, space (Catelan et al. 1998). 
With these assumptions, a typical cyclical term can be written 

£(ki,k 2 ) = A + Ai cos(6»i 2 ) + A 2 cos 2 (6 12 ) (10) 

where Q\ 2 denotes the angle between ki and k 2 , and 

!0 

A = — ci + c 2 

A 1 - ci ( — + — 

\ fc 2 fci 

A2 = |ci, (11) 
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and ci = I/61; ca = 62/61- Through these relations, the projected bispectrum will depend on the bias parameters b\ and 62- 
Using this, we will now calculate the theoretical expression for the projected bispectrum in spherical harmonics in the 
mildly nonlinear regime. With substitution Q we find 

I(n,r 2 ,r 3 ) = /i2 + I23 +I13 (12) 

Using the properties of spherical harmonics in Appendix A [equation ([i""]), the orthogonality relation, ( |4l| ) and (0)], we 
obtain that a typical cyclical term is: 



J 12 = (4tt)M tfki(fk 2 /C(ki,k 2 )P(fti)P(A: 2 ) X 
Jo 

^ i il i il (fcir 8 )y i *r(-a 1 )y i n I 1 (^) E ^ 2 ^ 2 (fc 2 r 3 )y i :; 2 (-^ 2 )y L " 2 2 (^3)- (13) 



We can now write 



dny™ 1 (fi)y™ 2 (fi)y™ 3 (n) = h™? , (14) 

4tt 



which can be expressed in term of Clebsch-Gordan coefficients, and is non-zero only if the following symmetry conditions are 
satisfied: 

• m 2 + m 3 = m\ 

• l\ + 62 + (.3 — even 

• satisfy the triangle rule 

From equation ([j~3|) we thus obtain: 
I dn ri dn r2 dn r3 Y™^n ri )Y™^n r2 )Y™ 3 (n r3 )h2 = (4tt) 4 f dk 1 dk2k 2 1 k 2 2 P(k 1 )P{k2)F(r 1 ,r2,r 3 ,k 1 ,k2) , (15) 

J 4-7T J 

where 

F(n,ra,rs,*i,fca) = / dn kl dSl k2 (A + A x cos8 12 + A 2 cos 2 i2 ) / 1+£2 (-^("i+^+^lj^ (fciri)^ 2 (fc 2 r 2 ) X 

E (fci-3)j, 7 (fc2r 3 )y; i - mi (n fcl )y/ 2 " m2 (n fc2 )y™ 6 (-" fcl )y/ 7 m7 (-^ 2 )^™r 7 ( 16 ) 

and F can be written as Fo + Fi + F2 where Fo involves the term Ao etc. 
The Fo term is easily calculated: 

Fo = A t 2( ' 1+ ' a) j tl (k iri )j e2 (k 2 r 2 )ji 1 (fcir 3 ) ji 2 (k2T3)Hi 3 m e lir 2 ( 17 ) 

and therefore satisfies the symmetry rules. 
For Fi and F 2 we exploit the fact that: 

COS 612 = Pi (COS 012 ), (18) 



2 

cos 



012 = - [2P 2 (cos0 12 ) +Po(cos0i 2 )] (19) 
and use the addition theorem for spherical harmonics: 

rt 

P4cos0 12 ) = ^ T E Y™(n kl )Y: m (n k2 ). (20) 

m — — n 

The Pi term then becomes: 
4tt 

T 



Pi = A^i^i-ir^j^rM^) e ^ + ' 7 ^(fci^)j, 7 (fc2fc 3 )H73 m 47r 7 ^r 6M < f 7v m7 - ( 2i ) 



It is easy to see that ma + m7 = ~m3. To demonstrate that the symmetry conditions are all satisfied, consider the following 
part of eq.(^l|): 

Enj-m 3 mg TT17 -Tj — tn\ —rriQ M. s%_j%fcm,% — rrif ( r ) r )\ 

mQinj M 
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let's introduce a new quantity h™ 1 ^ 2 ^ 3 that is symmetric for any permutation of the columns (J 1 *) - It is clear that: 

T_,m 1 m 2 m 3 / , \ mi . — mi m 2 m 3 fnn\ 



The quantity h™ 1 ^ 2 ^ 3 can be written in terms of the 3J symbols: 



KXT = V(2*i + l)(2£ 2 + l)(24 + l)/(47r) £ £ ) {^J 3 ) . (24) 
Equation ( J2l] ) therefore contains the following multiplicative term: 

V m 6 m 7 J \m\ — m e M / \— M m 2 -my / ^ x / V m l m 2 m 3 j \f 7 f e 1/ ! 



mg ni7 Af 

where ^ f } denotes the 6J symbol. In the last equality we used eq. (4.88) of Sobelman (1979). The properties of the 
3J symbol ensure that the Fi term satisfies the symmetry conditions. 
Similarly for the F 2 term we obtain 

F 2 = A 2 i^ +l2 {-l) m ^ m3 ^j tl {k 1 r l )j l2 {k 2 r 2 ) H)*"* 'jt 9 (kir 3 )j tT (k2r 3 ) x 



nj — m3m,Q m-j 

n £* £a £7 



9 

^ nj — m 1 M — rn q n/Mm.2 —m^ . ni — rn\0 — rn,Qnj0rn2 —mj 
K il 2 E G n 2 £ 2 £7 + £\ £q n £ 2 £7 



r ti Z £q Z t2 £7 «l U «6 U £2 ^7 

(26) 

The second term in the square brackets does not present any problem, in fact it is nonzero only if £\ = £e, £2 = £7, rn\ = m^, 
m 2 = 7717, and satisfies the symmetry conditions. Similar methods to those above complete the symmetry considerations. 

Factorising the Wigner 3J symbol, and collecting terms together, we find the expression for the angular bispectrum, as 
a sum of cyclical permutations: 

S W3 = B12 + B13 + B23 (27) 
where, writing ^t(k) = ~p J drr 2 je(kr)ip(r), 



1 16 / (24 + l)(% + 1) (24 + 1) f^^+t,^ 
n 3 ty y (47r) 



Bia = =s—\ (1 Yi I dk 1 dk2i Zl ^ 2 kik^P(k 1 )P(k 2 )^eAki)^t 2 (k2)x 



J2 ^ +£7 (-l) £ ^(fci,fc 2 )(% + 1)(2* T + l)p / drr^(r)j, 6 (fc 1 r)^ 7 (fc 2 r) «« „<) (f, 2 * <) («» «« {£ £ «» } (28) 
where Be(ki, k 2 ) for f = 0, 1, 2 are: 

B (k\, k 2 ) = ^-ci + c 2 

\k 2 kiJ 

Ba(fti,*a) = IfCi, (29) 

and the sum X)« 6 <? 7 exte nds over £ = 0, 1, 2; l % =t x -t £1 + £; £7 = £ 2 - £ £ 2 + £. 

The above expression can easily be generalized for any 3D bispectrum. In fact, since a) the bispectrum is non-zero only 
if the three k vectors form a triangle, b) the bispectrum does not depend on the spatial orientation of the triangle (isotropy) 
and c) a triangle is completely specified only by the magnitude of two sides and the angle between them, the bispectrum 
can always be expressed as a sum over three cyclical terms each involving only the modulus of two fe-vectors and the angle 
between them: 

fl(ki,k 2 ,k 3 ) = -F(fei,fe 2 ,0i2). (30) 
Each of the cyclical terms can therefore be expanded as: 

n 

.F(fci,fe2,0i2) = P(fei)P(fc2)^B J e(fei,A; 2 )P < (cos0i2) (31) 

1=0 

where now P(ki) is an arbitrary function of |k; |, the coefficients Be can depend on any combination of jki| and j k2 1 and the 
sum over £ should in principle go to infinity, but in practice will be truncated at n. 

We find that the exact expression for the projected bispectrum Be 1 e 2 e 3 is still given by equations (Effl and (jig) ) where 
now the sum over £ goes up to n. 

This, with equations ( 271) and (Eq), is the major new result of this paper. 
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2.1 Applications 

Equation (|2^) has therefore much wider applications than the second-order gravitationally induced bispectrum considered 
so far. The mathematics developed for this purpose can be straightforwardly applied to CMB and gravitational lensing 
studies. The gravitational fluctuations and cosmological structures along the path of the last-scattering surface photons, 
distort the CMB signal mainly through gravitational lensing, the integrated Sachs-Wolfe effect (Sachs & Wolfe 1967), the 
Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1980), and through the Rees-Sciama effect (Rees & Sciama 1968) and other 
second order effects (e.g. Mollerach & Matarrese 1997 and references therein). In particular, if the primordial fluctuations 
were Gaussian, many of these effects can introduce non gaussian features in the CMB signal. The bispectrum is a powerful 
tool for detecting these effects to probe the low-redshift Universe (Goldberg & Spergel 1998, Spergel & Goldberg 1998). 
Contributions to the CMB bispectrum induced by secondary anisotropies during reionization (Cooray & Hu 1999), non- linear 
gravitational evolution (Luo & Schramm 1994; Mollerach et al. 1995; Munshi, Souradeep & Starobinsky 1995) and foregrounds 
(e.g. Refregier, Spergel & Herbig 1998) imprint specific signatures on the CMB bispectrum which need to be subtracted from 
the signal in order to be able to test the gaussian nature of primordial fluctuations. 

On the other hand, primordial fluctuations can induce nonzero bispectrum in the CMB, that encloses information about 
the physical mechanism that generated them (e.g. Falk, Rangarajan & Srednicki 1993; Luo & Schramm 1994, Gangui et al. 
1994, Mollerach et al. 1995, Gangui & Mollerach 1996, Wang & Kamionkowski 1999, Gangui & Martin 1999). The evaluation 
of all these contributions to the observed CMB bispectrum requires calculation of an integral as in our equation (5) and (6), 
where r 2 ip{r) is replaced by an appropriate weight function. 

In the local Universe, gravitational lensing provides a direct probe of the mass fluctuations. The study of Fourier space 
correlation functions of the gravitational weak shear and convergence field is still in its infancy, but it is potentially fruitful: it 
could give us detailed knowledge of the correlation properties of the projected mass distribution (e.g. Bernardeau et al. 1997 
and references therein, Munshi 2000). 

In the present paper we will use equation ( p8| ) together with (^) as an exact expression for the second-order perturbation 
theory bispectrum of an angular catalogue with selection function ij)(r) and galaxy power spectrum P{k), assuming a local bias 
model with parameters b\ and fe- In principle, one can estimate the angular bispectrum from a projected galaxy catalogue, 
and use likelihood methods to constrain the bias parameters, which enter through the Be terms. In Section 3, we compute the 
likely errors from such a study, to determine if it is worthwhile to undertake such an analysis with current catalogues. Before 
we do so, it is worth noting that, in the current form, it is very expensive to compute: in the following subsection we rewrite 
it in a form more suitable for practical evaluation. 



2.2 Practical evaluation of Be 1 e 2 e 3 

From a computational point of view it is possible to speed up the calculations considerably (and consequently make the 
problem computationally manageable) by rewriting equation (^) in terms of the function Q\ defined as: 

ej,(*j,r) = J dk^ h {k)k 2 P{k)j t] {kr)k q (32) 

where q = —1,0, 1 and {i,j} = {1,6} or {2,7}. This function can be evaluated and tabulated in advance to speed up the 
analysis. With this definition, we can write the components of the angular bispectrum as 




x 



J drsrl^rs) (| Cl + c 2 ) 9?^, r 3 )e? 2 (^, r 3 ) ^ °) (§■ t o°) & 0* o°) {£ % 0} + 

cx £ [e7>,r 3 )e+>,r3) + e+ 1 (4,r 3 )e7 2 1 ^,r 3 )] t cO (? «f o l ) ^ ^ o f7 ) {% %^} + 

fl-K« 6 <«l + 1 
£ 2 ~l<e r <i 2 + l 

A C1 ®U^rs)el(£r,rs) (* * *) (* * „*) (fr *• *) {% % f) . (33) 

t 1 -2<e 6 <t 1 +2 
e 2 -2<i 7 <e 2 +2 

Note that this analysis is appropriate for all-sky coverage, and ignores shot noise. This is a good approximation for the high 
surface-density catalogues such as APM, in the range where perturbation theory is valid. Estimators for noisy data and partial 
sky coverage are presented in Heavens (1998) and Heavens (2000), see also Gangui & Martin (2000). Note also that numerical 
codes can run into difficulties when computing the spherical harmonic expansion and 3J symbols for high I. In Appendix 
C we give asymptotic expressions at high i for the 3 J symbols that are easily evaluated, and we present a way to calculate 
spherical harmonics fast and accurately at high I. 
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3 ERROR ANALYSIS FOR THE BIAS PARAMETER 

The spherical harmonic bispectrum in second-order perturbation theory is a known function of the galaxy power spectrum, and 
depends on the bias parameters b\ and 62 through So, Bi, B2 (equation^). Equation ( |28| ) relates therefore two measurable 
quantities (the spherical harmonic bispectrum of galaxies and the galaxy power spectrum) via the unknown bias parameters 
61 and 62. The 3D power spectrum may be obtained from the projected catalogue either by deconvolution of the angular 
correlation function or the angular power spectrum (e.g. Baugh & Efstathiou 1993; Baugh & Efstathiou 1994). In practice 
this is done in the small-angle approximation. In Appendix B we show that this is perfectly adequate for the power spectrum. 
We therefore have a full prescription for the angular bispectrum in terms of observable quantities and parameters which we 
wish to measure [equation ( ps| ) or the computationally manageable equation (p3|)]. The problem is therefore suitable for a 
likelihood analysis to extract the bias parameter. Such a programme is a major undertaking, so it makes sense to compute 
the expected error on the bias parameter first, to see whether the programme is likely to succeed. 



3.1 Likelihood analysis for ci and C2 

We assume we have full sky coverage unless otherwise stated. We define the quantity^ 4=Re[a^ 1 a^ 2 aJ™ 3 ] with £i, m, such 
that (mjmjmg 3 ) 7^ 0. For a given triplet £1,^2,^3 there are (2ti + 1)(2^2 + 1) distinct d a . From d a we can build the unbiased 



estimator of Bi lt2 i 3 , D a = „ e ° ^ 



Since D a is unbiased, any combination D a = , w a' /(^2 a , w a >), where w a i is a weight, is also an unbiased 
estimator of Bi 1 i 2 i 3 . 

The optimum weight w a that minimizes the variance {(D a — Bi 1 e 2 e 3 ) 2 ) is w a = l/o^ = („, m 2 m 3 ) /o"d Q (cf Gangui 
& Martin 2000). The minimum variance estimator is 

Z_/mi m2m3 — — 

D « = 77. ,'°,^2 (34) 

y 



mi m 2 m 3 j 
^ 



1 m 2 m 3 

The variance of d a does depend on m, but only weakly. There is a leading term, independent of m, proportional to 3 

ft, ■ £■ \ ft £ £ \ 

angular power spectra (C ei Ce 2 Ci 3 ), plus a sub-leading term proportional to B tiljlk Bi plq i T \ mi ,^ mfc J ^ q m *" r J, where 

{i, j, k,p, q, r} is a permutation of {1, 2, 3, 1, 2, 3}. If we ignore the m-dependence of this last term, then the estimator simplifies 
to 

D a = ^ d a (m 1 m 2 m 3 3 ) ! i = 1, 2, 3. (35) 
Strictly it is not the minimum variance estimator, but it is not far from it, is much simpler, and is unbiased. 



3.2 A priori error for the bias parameter 

Since the quantity (a^ 1 a^ 2 a^ 3 ) can be factorized as in equation (^) it is possible to evaluate the expected error on C\ 
estimation by approximating the variance by its leading term, neglecting shot noise and by considering uncorrelated data, 
obtaining: 

a <* { dcl>- 2^ Gll C i2 C l3 2^ N ti (m,) W 

where C denotes the likelihood function. 

The quantity Ne i (mi) denotes the number of terms like Ci 1 Ce 2 Ce 3 present in the covariance. It depends on the configu- 
ration i.e. on the choice of the triplets of £'s. It is useful to notice here that in the absence of Ni i (mi) we have 



E 



(37) 



Equation ( |36[ ) assumes that the covariance matrix is diagonal, this means that different bispectrum estimators (a™ 1 a^ 2 a™ 3 ) 
are uncorrelated. In the case of a survey with full sky coverage, similarly to the three-dimensional case treated in MVH97, 
the covariance matrix is well approximated by a diagonal matrix if each a™ appears in one estimator only. However, in the 
presence of a mask, different a™ are correlated, therefore (pfl) might no longer be valid. 



* Re [x] denotes the real part of the complex number x. 
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For an order-of-magnitude estimation of the expected error on the bias parameter, let us consider only equilateral 
configurations (i.e. configurations where £1 — £2 ~ £3), and assume full sky coverage for a survey with the APM selection 
function. It is easy to estimate the error achievable on ci using equation ( p^ ) and considering that second order perturbation 
theory should hold up to £ — 35. This choice is justified by the following argument: in the three-dimensional galaxy distribution, 
second-order perturbation theory breaks down at k ~ 0.6 (Mpc h^ 1 )^ 1 (cf MVH97, although it depends on the power spectrum 
slope and can be smaller, see e.g. Scoccimarro et al. (1998)) that corresponds to a scale of the order of 10 Mpc hT 1 . At the 
medium depth of the APM survey (335 Mpc /i" 1 ), this subtends an angle of about 0.03 radians (in agreement with the findings 
of Gaztanaga & Bernardeau 1998), corresponding with I ~ 33. This order of magnitude calculation yields an estimate for the 
error on Ci of about ±3.5, which is not really encouraging. However, this is only an order-of magnitude calculation: a more 
rigorous treatment is implemented in the next section. 



3.3 The choice of the triplets 

As already discussed in MVH97, the choice of the triplets to evaluate the bispectrum is very wide, but, speed and memory 
considerations force one to simplify the analysis by ensuring that the covariance matrix is diagonal, for a full sky survey, this 
can be achieved by ensuring that each I appears only in one triplet. The choice of the ratio between the ts (the shape) of a 
triplet, is influenced by the behavior of the bispectrum: triplets with the same shape give an almost degenerate information 
on ci and C2, in practice each shape can constrain a linear combination of ci and C2: the likelihood will be aligned along 
a straight line in the ci,C2 plane. The best choice to try to lift this additional degeneracy is to combine the likelihood for 
equilateral triplets (£1 — £2 = £3) with the likelihood for degenerate triplets {£\ = £2 and £3 = 2£\). 



3.4 Covariance 

To perform a likelihood analysis we need an expression for the covariance matrix for our estimator Bi 1 i 2 i 3 . It is easy to verify 
that, if each £ appears only in one Bi 1 i 2 i 3 then the B[ 1 i 2 i 3 are uncorrelated, that is: 

{B ei i 2 i 3 Be 4 e 5 e 6 ) = (-Bfi^faX- 6 ^^) (38) 
if £i 7^ £j, where i = 1, 2, 3 and j = 4, 5, 6. This means that the off-diagonal terms of the covariance matrix are zero. In fact: 

{B ei i 2 i 3 B ei e 5 e 6 ) = ^ ^ {a^a^a^a^a^a^ 6 } m 2 2 J 3 ) m J 5 m 6 J (39) 

m\m 2 Ta 3 771477157716 

where we used the fact that a™ = (— l) m aj m . Analogously to MVH97, the quantity {a™^ aj^ 2 a™ 3 a™^ a™ 5 a™ 6 6 ) can be split 
into: 

a) 15 cyclical permutations of the kind (a^ 1 a^ 2 )(a^ 3 a^ 4 )(a^ 5 a^ 6 ) that are all zero if £i 7^ £j, where i — 1,2,3 and 
j = 4, 5, 6, 

b) 1 term 

Be 1 i 2 e 3 Bi 4es £ 6 ^ ^ (m, m 2 ^3) (m 4 m 5 m 6 ) = (Be 1 t 2 e 3 } {Be 4 t s i 6 } 

771177797773 m^m^mQ 

c) 9 cyclical permutations of the kind: 

d d (** l 'i £ " \ (t'i l i l ' k \ (ii 12 t 3 \ (tt e 5 e 6 \ 

"titjth / j / j \mi m'. m k J ym'. mj m' k J V m l m 2 >™3 / V m 4 »™5 m 6 ) 

where k is any permutation of 1,2,3 and i'j'k' denotes any permutation of 4,5,6. These terms in c) are all zero unless there 
are repeated £ in two different D a . 

The term in b) cancels when subtracting the mean to obtain the covariance. Let us now consider the diagonal terms of the 
covariance matrix: these are given by equation (^) with the following identities for the indices: 1 = 4, 2 = 5, 3 = 6. For 
symmetry considerations we can restrict ourselves to consider £\ < £2 < £3. In the case where £\ < £2 < £3 we have: 

in a) only one term surviving, giving Ct 1 Ce 2 Ct 3 ■ 
in c) using (||), SBj^^ 

In the case where l\ = £2 < £3 we have: 

in a) Cf 1 Ci 3 + ^2 mm , (m-mo' 3 ) (m' ^m'^o)] m ^ ne particular case where £3 = 2£\ (degenerate configurations) as 
shown in Appendix C (eq. ^32j) we obtain to very good approximation: Cj 1 Ci 3 2 + yj (2tt£i)/(1 + 4£i) 
in c) 5-B? , 1, r „ 
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Figure 1. The power spectrum in spherical harmonics for the simulated projected catalogue. The points are the Ci measured for the 
catalogue, the solid line is the underlying power spectrum, and the dashed line is the underlying linear power spectrum. The underlying 
power spectrum is obtained from the 3D one by applying the selection function and the full-treatment projection. Deviations from linear 
theory are already evident at I ~ 40. 

For equilateral configurations where l\ = £2 — £3 — I we have: 

in a) Cf [6 + 9£ mm , ( e m J m e ) (t'-m^o)] which, to a very good approximation, is ~ Cf [6 + 9 X 1.15/(2^ + 1)] 
in c) 9Bj u 

3.5 Likelihood analysis of a simulated catalogue 

We created an all-sky catalogue with the APM selection function (f>(r) oc r _0 1 exp [(-r/335) 2 ] (e.g. Maddox, Efstathiou & 
Sutherland (1995), Gaztanaga & Bernardeau (1998)), by replicating an N-body simulation of 500 Mpc h~ side and then 
sparsely sampling and projecting to the plane of the sky 2303636 particles (galaxies for our purposes) accordingly to the APM 
selection function. The simulation was supplied by J. Peacock using the AP 3 M code (Couchman 1991). The characteristics 
are: 128 3 particles, CDM transfer function (Bardeen et al. 1986), P = 0.25, Q. = 0.3, A = 0.7, evolved to ag=l to best fit 
the observed cluster abundance; this choice gives also a good fit to the COBE 4-year data (e.g. Tegmark 1996). The rate of 
sampling used ensures that the probability of selecting the same particle in the simulation more than once in the replicated 
box is negligible. 

Figure |l| shows the angular power spectrum for the simulated projected catalogue. The solid line is the underlying power 
spectrum obtained from the 3D one by applying convolution with the selection function and full-treatment projection; the 
dashed line is the underlying linear power spectrum. Deviations for linear theory are already evident at £ ~ 40. In order to 
lift the degeneracy between ci and C2 we consider equilateral and degenerate configurations. The likelihood for equilateral 
configurations is shown in Figure ^; it does not give a strong constraint on the bias: perturbation theory breaks down at 
I ~ 40: up to I < 40 there are only 18 independent equilateral triplets^. The likelihood for degenerate configurations gives a 
better constraint on the bias. Perturbation theory for this configuration breaks down where the short I is £ = 40. Likelihood 
contours for degenerate triplets configurations where the short £ is 20 < I < 40, are shown in Figure |j| Since the likelihood 
for equilateral configurations is quite broad, even combining it to the likelihood for degenerate configurations does not modify 
Figure ^| sensibly. 

From this analysis we can conclude that from a two-dimensional galaxy survey with the APM selection function, even if it 
is an all-sky survey, it is only possible to constrain a combination of the linear and quadratic bias parameter, or, alternatively, 
if we assume (without justification) that the bias is linear (i.e. &2 = 0), 0.7 < ci < 1.4 or 0.7 < 61 < 1.4, at 68% confidence. 




t 18 is the number of independent D a as defined in eq. 35 with i$ = I < 40, i = 1,2,3. The presence of the 3-J symbol for equilateral 
configurations requires I to be even, £ = 2 is discarded because it would be contaminated by the galaxy quadrupole. 
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Figure 2. Likelihood contours for equilateral triplets configurations. The two levels (where the downhill direction is indicated for clarity) 
are the 1—cr and 3-cr confidence levels and the + indicates where the true value for the parameters lies. Perturbation theory breaks down 
at £ ~ 40 here £ up to 42 are considered even though the likelihood only for the triplets between £ = 40 and £ = 42 includes the true 
value just within the boundary of the 3 — cr level. This configuration does not place strong constraints on the bias parameter. 




Figure 3. Likelihood contours for degenerate triplets configurations. The two levels are the 1— a and 3-cr confidence levels and the + 
indicates where the true value for the parameters lies. Perturbation theory breaks down at ^ s hort 

= 40. 



4 CONCLUSIONS 

In this paper, we have presented the formalism for translating the 3D bispectrum of a sample population into the angular 
bispectrum in spherical harmonics. As discussed in section 2.1, this method has applications in a variety of areas, such as 
microwave background studies, gravitational lensing and analysis of angular galaxy catalogues. We have investigated the last 
of these in detail in this paper: since the bispectrum is a measurable quantity, and its theoretical expression depends on 
measurable quantities via the unknown bias parameter, it is possible to extract the bias parameter via a likelihood analysis. 
We have therefore investigated its use as a tool for measuring the bias parameter for projected galaxy surveys. In principle, 
it is an alternative method to using 3D galaxy redshift surveys, without the complicating effects of redshift distortions and 
higher shot noise. Recently, other methods based on second-order perturbation theory have been proposed to measure the bias 
parameter from 2D galaxy catalogues (see e.g. Frieman & Gaztanaga 1999; Fry & Thomas 1999). Frieman & Gaztanaga (1999) 
studied the reduced 3-point correlation function on the sky. The error analysis in real space is more complicated because of 
strong correlations between the estimates. Frieman & Gaztanaga conclude that b\ <C 1.5 or so, giving a comparable error to 
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our analysis (section 3.5) if 62 is assumed to be zero. We emphasize that allowing a non-zero quadratic bias term opens up a 
wide range of acceptable linear bias parameters. 

The analysis of Fry & Thomas (1999) is closest to ours. They consider the bispectrum, but present results in the small- 
angle approximation only. They do go some way in writing down the general expression for the angular bispectrum in spherical 
harmonics in terms of the 3D bispectrum. In this paper, by expanding the (general) dependence of the 3D bispectrum on angle 
between wavevectors in Legendre polynomials, we were able to derive a practical general relationship which is computable 
with few numerical integrations. 

We have calculated the expected error on the linear bias parameter from an all-sky catalogue with a selection function 
similar to the APM survey. We find that the results are not encouraging for projected catalogues, and that it is preferable to 
undertake a bispectrum study of 3D galaxy redshift surveys such as the AAT 2dF or the Sloan Digital Sky Survey, using the 
methods discussed in MVH97 and Verde et al. (1998). Tests on simulated projected catalogues confirm our analytic findings. 
In a similar way as for 3D surveys (see discussion in MVH97), one can reduce the errors by subdividing the sky [J but by 
modest factors which do not change this basic conclusion, and at the expense of considerable complication in the analysis 
through window convolutions. 
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APPENDIX A: SPHERICAL HARMONICS 

Spherical harmonics are the natural basis for describing a two-dimensional random field on the sky. The definition we use 
here is expressed in terms of the associate Legendre functions Pp (cos 9): 



— vmrn ,\ / (2^ + l)(t ~ m)l „ m , „, m 6 f (-i) m for m > /,„% 

1 (") = 1 = y 4^ + m)! Pl {cOSd)e X b for m < ~ ( 40 ) 

where £ and m are integers and £ > 0, —I < m < £. Their orthogonality relation is: 

dQY e ^ (n)YiT 2 = 8? lta 5* ima . (41) 
Any two dimensional pattern f(Q) on the surface of a sphere can be expanded as: 
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f(n) = J2°-TYr*{fy (42) 

where 

af = J ' dQY e m (n)f(Q). (43) 

Useful relations involving the spherical harmonics are: 
y f m *(fi) = (-l) m Yf m (a) , Y e m (-0) = {-l) m+e Y e m (Sl). (44) 
and the identities: 

^y/ m (fi)Y/"(fi') = 6(n - n') (45) 

exp(ik • r) = 4iY^2i e j e (kr)Y e * m (n k )Y e m (V r ) (46) 

lm 

e 

Pt(fx*B) = ^^WW (47) 

m=-t 

where 8 denotes the angle between the vectors r and k. The latter is the addition theorem for spherical harmonics. 

APPENDIX B: THE POWER SPECTRUM AND THE SMALL-ANGLE APPROXIMATION 

The power spectrum of a 2-D distribution on the plane on the sky is given by the set of Ce denned as: 

{aTap *) = Ci6ui6* m i. (48) 

The corresponding angular two-point correlation function can be expanded as 

C ( 9 ) = T- V(2«+l)CeP/(cas0) (49) 

4-7T L — ' 

e 

with inverse relation is: 

Cn = 2n J C(8)Pe(cos6)dcos(6) (50) 

On the other hand, for a plane two-dimensional distribution with power spectrum P2d(k) the two-point correlation function 
is: 

"OO />27T 

2 . 1 / I „ , S , „ , X , , , I 



w(8) — I P2d{k) exp(in ■ 6)d a = / / P2d(k) cos(n8 cos cj>)d4>KdK — j-—. - / P2d(k) Jo(K,8)Kdn (51) 

( 27r ) J ( 27r ) Jo Jo ( 2n > Jo 



In the small angle approximation, i.e. in equation ( |49[ ) for small 8, we have that P^(cos#) ~ Jo[(l + 1/2)6], but, since small 
angular patches restrict us to high I, Pt(cos8) ~ Jo(£8). We can therefore conclude that in the small-angle approximation 

t — >k; C e — yP2o{l) (52) 

For a real angular catalogue the two-dimensional galaxy density in the sky is obtained as follow. Let the true three-dimensional 
galaxy density field be p(r) and the selection function be ip{r), normalized here to J drr 2 ip(r) = 1. It is straightforward to 
obtain an expression for the angular power spectrum given the three-dimensional one: 

{a mi a m2 *} = ( W^I dkk2p ^ [J drr2 ^( r )jii( k i r )] 2 if m i = ~ m 2 and h = t 2 ^ 
ei t2 y otherwise 

In the small-angle approximation this is (Kaiser 1992, Buchalter, Kamionkowski & Jaffe 1999): 

P*d(k)= drP 3D (K/r)ip 2 (r)r 2 (54) 
Jo 

Also in the presence of the selection function we can check that mapping ( p32"| ) is valid and we can asses the limit of validity 
for the small angle approximation for the power spectrum. 

Assuming the APM selection function 4>(r) oc r~° ' 1 exp[(-r/335) 2 ], and a CDM power spectrum (Efstathiou, Bond & 
White 1992) with V — 0.25, we compared the angular power spectrum obtained with the exact projection as in equation 
( [si]) and with the small angle approximation [equation Q . The result is shown in Figure [| the small angle approximation 
introduces an error smaller than 3% for £ > 20. 
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Figure 4. The small-angle approximation for the angular power spectrum works very well for I > 20 introducing an error smaller that 
3%. 

APPENDIX C: USEFUL FORMULAE FOR THE HIGH £ REGIME. 

Library routines dealing with spherical harmonics at high £ can sometimes fail. In this appendix, we present some asymptotic 
results which can avoid problems. 

The 3J symbol (£ lf 2 e §) may be written as: 



f<iWt_ t.sL l (L-2h)\(L-2£ 2 )\(L-2e a y. (L/2)! 

^° 00 ^ 1 j V + (L/2-£ 1 y.(L/2-e 2 y.(L/2-e 3 y. (> 



(56) 



where L = l 1 +l 2 + £3. 

In the special case where £\ — £2 — £ and £3 = 21 this becomes: 

2[r(2i)] 2 

VViT4Z[r(£)] 2 yr(47) 

Numerical routines to calculate the 3J symbols usually encounter problems for large £. We therefore evaluated an approx- 
imation based on the Stirling approximation: i.e. n\ — T(n + 1) and (Gradshteyn & Ryzhik 1965) 

r(z) ~ e- z z z - 1/2 (2n) 1/2 for large z. (57) 

This approximation for the F function is quite good, in fact it introduces an error of only 4% at z — 2. Using this approximation 
we obtain for eq. (fiB]): 

/ i\L/2 9 l/4 r~ tL/2+1/2 

{0 > v/2i (X + l)V2+i/2+i/4 p/2 _ fl )(L/ 2 - fe)(L/2 - f 3 )]i/i ■ 1 ; 

This approximation introduces a small error of a few percent at £ ~ 20. For equation (|5(j) we obtain: 

1 / A \ 0.25 

When the m's are non-zero it is still possible to find a simple expression for the 3J symbol in special cases. For example if 
£3 = £1 + £2 we have: 



/*1 «2 
\m1m2 —mi —77X2 



e 1+ i 2 \ _ f i\li-<2+mi+™a / (2li)!(2l 2 )!(li + £2 + mi + m 2 )!(li + £2 - mi - m 2 )T~ 

™ 1 - m2 ' ,_l J V (24 + 2< 2 + 1)!(4 + mi)!(< 2 + m 2 )!(t - mi)!(fe - m 2 )! ( 1 

In the special case where £\ — £2 = £ the previous expression can be further simplified and approximated -using (^)- by: 
(it 2i \ , n m 1+ m 3 (12tt) 1/4 / (2f + m 1 +m 2 )!(2f-m 1 -m 2 )! 

(mi77i2-mi-77i 2j / > 2MV3T+TV (i + mx)!(^ + m a )!(/-mi)!(/-m a )! 1 j 

In the calculation of the covariance matrix for "degenerate" configurations for the quantities B«2^ we came across with the 
following sum over m of a product of two 3J-symbols. An useful expression for it is the following: 
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\mi -mi ) \ 



\ 2l\ 
m 2 — m 2 J 



= 2 4i £[(2£- l)!] 2 

(i + «)r(4£)r(2£+ 1) 2 ~ (1 + 41) (4£-i)V 



m l > m 2 

For large £ using the above approximation for the Gamma function we obtain: 
(l + 4Qr(4Q (1+40 ( g ] 

this approximation works very well also at low ^, in fact introduces an error below 1% at £ = 6. 
The orthogonality relations for 3J-symbols are also widely used: 



(tl <2 Id. y _ I 
^ \ml m2 m3 / 

mi m2 1T13 



s 



and 



E 



1 ml ml m3 



) (™i 



1 



(2£ 3 + 1) 



(62) 



(63) 



(64) 



(65) 



When calculating the spherical harmonic coefficients for a galaxy distribution on the celestial sphere, numerical problems 
arise with the associate Legendre polynomials at high £. Routines based on the recurrence relations used for example by the 
numerical recipes routines (Press et al. 1992) fails at £ ~ 35 if m ~ £, at slightly higher £ for m <C Z. The asymptotic expansion 
for the associate Legendre polynomial (e.g.Gradshteyn & Ryzhik 1965): 

2 r(^ + m + l)cos[(^ + 1/2)61 -7r/4 + m7r/2] 



P™(cos ( 



r(^ + 3/2) 



+ 0(l/£) 



(66) 



V2 sin 9 

is valid for | £ |>| m |, 1 1 \> 1 and e < 6 < tv — e. 

The use of the recursive relations involve the calculation of (2£ — 1)!! that at high £ can create numerical problems. Using 
the following expression for the F function: 

r(n + l/2) = ^(2n-l)!! (67) 
and (^?j) for the T function for big argument we obtain: 

This approximation introduces an error of only a few percent for n ~ 10. 

When calculating the spherical harmonics at higher £, problems arise not only with the associated Legendre polynomials, 
but also with the part that involves the ratio of two factorials. A better way to calculate the spherical harmonics, fast and 
accurate to high £ is based on the algorithm proposed by (Muciaccia, Natali & Vittorio 1997). In essence the numerical 
problems can be avoided by defining the normalized associated Legendre polynomials XT: 



XT (cos 0) = 



21 + 1 (i - m) 
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The recurrence relation for AT is 
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with expressions for the starting values: 
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Am+i{x) = xV2m + 3XT l (x) 
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(69) 



(70) 



(71) 



(72) 



numerical evaluation can be greatly speeded up by noticing that the factor in ( frf] ) depends only on m and can therefore be 
calculated and/or tabulated for m < £ ma x only once. 
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